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ABSTRACT 

The decimetric luminosity of many galaxies appears to be dominated by 
synchrotron emission excited by supernova explosions. Simple models suggest 
that the luminosity is directly proportional to the rate of supernova explosions 
of massive stars averaged over the past ^ 3 x 10^ yr. The proportionality may 
be used together with models of the evolving 1.4 GHz luminosity function to 
estimate the global star formation rate density in the era z ^ 1. The local 
value is estimated to be 0.026 Mq yr“^ Mpc“^, some 50% larger than the value 
inferred from the Ha luminosity density. The value at 2 ; ~ 1 is found to be 
0.30 Mq yr“^ Mpc“^. The 10-fold increase in star formation rate density is 
consistent with the increase inferred from mm-wave, far-infrared, ultra-violet 
and Ha observations. 


Subject headings: cosmology: observations - galaxies: evolution - radio 
continuum: galaxies 


1. INTRODUCTION 

Two different mechanisms are believed to energize the relativistic electrons responsible 
for the decimetric synchrotron emission from galaxies. In classical radio galaxies and related 
objects, the energy source is ultimately a non-stellar central engine, perhaps a massive 
black hole. In normal disk galaxies and related objects, the energy source is ultimately 
stellar, perhaps the totality of supernova explosions over the past 3 x 10^ yr. Although 
the relative importance of the two mechanisms in some galaxies is debated, it is usually 
not difficult to determine the dominant source in any galaxy from observations of its radio 
morphology, and the brightness, morphology and spectrum of any optical counterpart. 

For galaxies in which the radio emission is energized by stars (hereinafter star-forming 
galaxies) it should be possible to infer the underlying rate of star formation from the 
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observed radio luminosity. Because the radio luminosity function of star-forming galaxies is 
well determined locally, and measurable out to ~ 1 , and because radio emission usually 
does not suffer absorption, the derivation of star formation rates from radio luminosities 
offers an important addition to other methods used to determine global star formation 
rates (e.g., Gallego et al. 1995; Rowan-Robinson et al. 1997; Pettini et al. 1997; Tresse & 
Maddox 1998; Treyer et al. 1998; Blain et al. 1998). On the other hand, the quantitative 
link between 1.4 GHz luminosity and the rate of star formation is rather tenuous, and 
further work is clearly needed to strengthen the astrophysical basis of the method. 

This note derives the star formation rate density for z ^ 1 from the observed radio 
luminosity function, and compares the result with other determinations of this important 
quantity. A Hubble constant Hq = 50 km s“^ Mpc“^ has been adopted, and previously 
published quantities dependent on Hq have been scaled if necessary. 


2. THE STAR FORMATION RATE AT 0 


Gondon (1984a, 1984b, 1988, 1989) has determined the local radio luminosity function 
for the combined population of spiral and irregular galaxies in the form 
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Here, L is the luminosity (W Hz“^) at 1.4 GHz, 0^ = ^ogf.{m)L(f){L) is the luminosity 
function per unit interval in log^(L) and per unit volume (Mpc“^), and Lq = 10^^'^"^ is the 
factor converting between SI units and astronomical units. The corresponding 1.4 GHz 
luminosity density, £ 1 , 4 , is given by 
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£14 = / L0(L)dL = 1.8 X 10^^ W Hz"^ Mpc'l 

J O 

The uncertainty in this estimate (Gondon 1989, Table 3) is a factor of 10^°'^. 


Gondon and Yin (1990) and Gondon (1992) have determined a relationship between 
the 1.4 GHz luminosity of a galaxy and its current star formation rate (SFR) which can be 
expressed in the form 


^FRi.4 = 
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The conversion factor in Equation (2) corresponds to an adopted initial mass function 
^/(M) ~ M“2-35 -with 0.1 < M < 100 M©. Equations (1) and (2) imply that the local value 
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of the global star formation rate density is 0.026 Mq yr“^ Mpc“^, with an uncertainty 
of the order of a factor of 10^°'^. The most luminous galaxy encountered in constructing 
the local luminosity function of the star-forming population has a 1.4 GHz luminosity of 
^ 10^^'^ W Hz“^ (Condon 1989, Table 3); its corresponding inferred star formation rate is 
just over 100 Mq yr“h 

3. THE STAR FORMATION RATE FOR ^ 1 

Sub-mJy radio surveys not suffering from confusion were hrst made by aperture 
synthesis telescopes in the early 1980s. The shape of the 1.4 GHz source-count distribution 
revealed by these studies provided compelling evidence for signihcant evolution of the star¬ 
forming galaxy population. As reviewed by Condon (1988), the source count distribution 
normalised to the Euclidean slope [i.e., the curve N{S)S~‘^'^] flattens for S' 1.4 ^ 5 mJy. 
This flattening cannot be interpreted in terms of the population of AGN-powered galaxies 
which dominate the source counts above 10 mJy, nor can it arise from a non-evolving 
population of star-forming galaxies having the local luminosity function presented above 
(Danese et al. 1987). The faint counts can be reconciled with the local luminosity function 
of star-forming galaxies by allowing signihcant evolution of the star-forming population 
in the interval z ^ 1. Danese et al. (1987; see also Rowan-Robinson et al. 1993; Hopkins 
et al. 1998) have shown that consistency between the local luminosity function and the 
faint source counts for the star-forming population can be obtained by a pure luminosity 
evolution model satisfying Pia{z) oc (1 + z)^. 

In studies of the entire radio source population, Condon (1984a, 1984b, 1988, 1989) has 
established consistency between the evolving luminosity functions of both star-forming and 
AGN galaxies and the radio-frequency source counts over six orders of magnitude of hux 
density. His model uses a density evolution factor g{z) and a luminosity evolution factor 
f{z) in the form 

0 m(L, z) = g{z)4>m[L/f{z), z = 0]. 

This is the ‘co-moving’ luminosity function; Condon adopted a cosmological model with 
Do = 1 and h = 0.5. Values of f{z) and g{z) taken from Condon (1984b, Table 1), combined 
with Equations (1) and (2), lead to the 2 :—dependence of the global star formation rate 
density illustrated in Figure 1. The curve is consistent with that which would be derived 
using pure luminosity evolution P{z) oc (1 -f zY. 

Although Condon’s (1984b) model establishes consistency between the evolving 
luminosity function and the source counts, it is not unique. Nevertheless, source-count 
models involving signihcant evolution of the population of sub-mJy galaxies hnd growing 
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support as studies of their optical counterparts show signs of vigorous star-formation. For 
example, Kron, Koo & Windhorst (1985) found that the optical counterparts of the sub-mJy 
population include a blue component, dubbed ‘B radio galaxies’, with peculiar optical 
morphologies suggestive of interactions and mergers. Benn et al. (1993) have shown that 
many faint radio sources appear to be more-distant counterparts of the IRAS star-forming 
galaxy population. This population probably overlaps with the ‘B radio galaxies’. Hubble 
Space Telescope images suggest that the majority of sub-mJy radio sources are starburst or 
post-starburst objects, embedded in luminous disk galaxies and triggered by interactions 
(Windhorst et al. 1995). 


4. COMPARISON WITH OTHER ESTIMATES 

Gallego et al. (1995) have estimated the global star formation rate at z ^ 0 from 
their observations of the local Ho; luminosity density. Their derived star formation rate 
density is 0.013 M© yr“^ Mpc“^, one-half the value derived above. While a factor of two 
uncertainty in the global star formation density derived by such different methods is not 
disturbingly large, it is of interest to explore the origin of the difference. Accordingly, Figure 
2 compares the 1.4 GHz and Ha luminosity functions and the corresponding distributions 
of the inferred star formation rate densities, both expressed on a common log{SFR) scale 
to facilitate comparison. 

Over the range where the Ha luminosity function is well determined 
(4 X 10^^ < Lffa < 4: X 10^® W or 0.4 < SFR^ < 40 Mq yr“^) the two luminosity 
functions cross, with the Ha distribution rising somewhat more steeply to lower luminosities 
and dropping slightly below the 1.4 GHz function at the highest observed point. At faint 
luminosities the extrapolation of the Ha luminosity function rises above the 1.4 GHz 
function, but the difference never exceeds a factor of 5 even at the faintest luminosities. The 
star formation rate density inferred from Ha is approximately twice the density inferred 
from 1.4 GHz observations, for values of SFR below about 40 Mq yr“^. 

The excess of the 1.4 GHz value of SFR relative to the Ha value therefore arises at 
high values of SFR. This is clear from Figure 2, which shows that the Ha luminosity 
function is Ri 10“® dex“^ Mpc“^ at SFR^ = 100 Mq yr“^ while at the same SFRi ^, the 
1.4 GHz luminosity function is a factor of 25 greater. Luminosity weighting amplihes this 
difference in the plot of the star formation rate density. As a result, the upper quartile of 
star formation rate densities lies above 25 Mq yr“^ according to the 1.4 GHz data, and 
above only 15 Mq yr“^ for the Ha data. Even more extreme is the high star formation rate 
tail which extends to values > 100 Mq yr“^ in the 1.4 GHz data, but is absent in Ha (see 
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also Cram et al. 1998; objects in this category include Arp 220 and NGC 6240). The excess 
SFR density in the 1.4 GHz data at high values of SFR more than outweighs the excess 
Ha SFR density at low values of SFR, and accounts for the somewhat higher global rate 
deduced from the 1.4 GHz data. 

Other estimates of the global star formation rate at low redshift {z ^ 0.3) have 
been made by Treyer et al. (1998) and Tresse & Maddox (1998). The former, based on 
balloon-borne observations at 200 nm, yields 0.006 Mq yr“^ (adjusted to h = 0.5) at 

0.15, and the latter, based on Ha spectra from the /-selected Canada-France Redshift 
Survey, yields 0.020 Mq yr“^ at z ze 0.2. The Ha value is compatible with that of Gallego 
et al. (1995), while the far-UV value is smaller, perhaps signihcantly so. Both are smaller 
than the 1.4 GHz value. 

There have been several determinations of the evolution of the SFR density beyond 
2 : 0.25 (e.g., Lilly et al. 1996; Madau et al. 1996; Connolly et al. 1997; Rowan-Robinson 

et al. 1997). All agree that there is a signihcant rise between z = 0 and z = 1, of a factor of 
order 10. The 1.4 GHz results are consistent with this trend. 

It should be noted that dust extinction leads to uncertainty in estimates of the star 
formation rate from UV and optical indicators. For example, Pettini et al. (1997, their Fig 
7) suggest that an upward revision of a factor of approximately 2-6 should be applied 
to uncorrected estimates of star formation rates based on ultraviolet observations. Their 
revised values he close to, but below, those estimated in this paper in the range 0.2 2 : ^ 1. 

Blain et al. (1998, their Figure 8) have estimated the evolution of the global star formation 
rate using sub-mm observations which, like radio observations, should be unaffected by 
extinction. Their estimates agree, to within a factor of two over the range z ^1, with those 
derived here. 


5. DISCUSSION AND SUMMARY 

In addition to the consistency with sub-mm observations noted above, there is close 
agreement between the evolution of the SFR density derived from 1.4 GHz data, and 
from far-infrared observations of the Hubble Deep Field by the Infrared Space Observatory 
(Rowan-Robinson et al. 1997). This result could have been anticipated, in view of the 
strong correlation between 1.4 GHz and far-infrared luminosities in local star-forming 
galaxies, and the viability of theories of the correlation which relate both luminosities to the 
number of massive stars in the galaxies (e.g. Gondon 1992). Estimates of the SFR density 
inferred from far-UV or Ha data (uncorrected for extinction) appear to be systematically 
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lower, by a factor of order 2 — 4. The discrepancy could be a consequence of a moderate 
amount of optical/far UV extinction, an interpretation that finds some support in studies 
of nearby galaxies which have revealed evidence of widespread extinction, including high 
attenuation in some objects (e.g.. Young, Kleinmann & Allen 1988; see also Pettini et 
al. 1997). However, it should also be noted that the ratio of the instantaneous luminosity 
of various SFR indicators to the short-term average of the SFR itself is a complex function 
of time (e.g.. Ho; is weighted strongly to the most massive stars’ lifetimes, while 1.4 GHz 
radiation is weighted to the more abundant stars lying just above the lower mass limit for 
supernovae). It is likely that factors other than extinction will also prove to be relevant in 
explaining the discrepancy between various SFR indicators. 

Further studies of the extinction in galaxies having high values of S'F’Ri .4 and low 
values of SFR^ could help resolve this question. Of particular interest are galaxies having 
SFRia ^ 25 M© yr“^. In such galaxies, the ‘shortfall’ in Ha emission is particularly large, 
and this might lead to obvious characteristic differences between the morphology in 1.4 GHz 
and Ha images. 

In Gondon’s (1984b) model the evolution of the star-forming galaxy luminosity function 
is constrained almost exclusively by the fit to sub-mJy 1.4 GHz source counts. This is 
not a particularly tight constraint, and as more redshifts become available for sub-mJy 
star-forming galaxies the evolution will become better dehned. Already, the studies of Benn 
et al. (1993), Rowan-Robinson et al. (1993) and Hopkins (1998) suggest that consistency 
arises between the local 1.4 GHz luminosity function of star-forming galaxies, the sub-mJy 
source counts and the available n{z) distributions if there is luminosity evolution according 
to (1 -f zY. This picture is consistent with results obtained at other wavelengths (e.g. 
Hacking, Gondon & Houck 1987). Since the 1.4 GHz data suggest that approximately 50% 
of star formation at z = 1 takes place in star forming galaxies with L 1.4 10^^ W Hz“^, and 

since such galaxies are not difficult to detect at radio frequencies (S'1.4 ~ 250/iJy), optical 
spectroscopy of complete samples of sub-mJy 1.4 GHz sources could make an important 
contribution to our understanding of cosmic star formation history in the era z ^1. 

It is worth noting that studies of the cosmic evolution of powerful radio galaxies (i.e., 
AGN-powered objects with Pi. 4 ghz > 10^® W Hz“^) suggest that the population undergoes 
pure luminosity evolution according to ^ 1 . 4 ( 2 :) oc ( 1 - 1 - 2 :)^ out to z ^ 2, with a luminosity 
dependent redshift cut-off at higher 2 : (e.g. Dunlop & Peacock 1990). As emphasised by 
Dunlop (1998), these studies provide a robust determination of the evolution of the 1.4 GHz 
luminosity density, despite remaining uncertainties in the high-;^ evolution of the luminosity 
function itself. It is a remarkable fact that the luminosity evolution of the AGN-powered 
population appears to be quite similar that of the star-forming population (Danese et 
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al. 1987; Dunlop 1998). Dunlop (1998) advances the suggestion that the star formation 
rate and the AGN (black hole) fueling rate might both be related to the global rate of 
gravitational accretion/condensation. 

The use of 1.4 GHz radiation to determine star formation rates offers some advantages 
over alternative methods. For example, 1.4 GHz radiation does not usually suffer absorption, 
and the local 1.4 GHz luminosity function of star-forming galaxies is well determined over 
a wide luminosity range. Gompared with mm-wave and far-infrared (FIR) indicators, 
which share these advantages over optical and UV indicators, deep 1.4 GHz observations 
have superior sensitivity, and offer astrometric precision sufficient to make direct optical 
identihcations of even the faintest sources. The main weakness is the tenuous linkage 
between 1.4 GHz radiation and the astrophysics of star formation. Spatially resolved 
observations probing the relationship between decimetric radio emission and star forming 
processes in selected galaxies (e.g. Duric et al. 1995; Marsh & Helou 1998; Smith et 
al. 1998) tend to support the model used above, but further work is clearly needed to 
strengthen the astrophysical basis of the method. 
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Fig. 2.— Luminosity distributions (left panel) and star formation rate densities 
(right panel). Symbols denote the indicator used, and the redshift of the 
determination. The dash-dot curve depicts extrapolations of the local Ha and 
1.4 GHz distributions beyond the brightest and faintest detected objects. The 
abscissae, labelled “star formation rate”, represent a common scale for 1.4 GHz 
and Ha luminosities, corresponding to Li 4/2.35 x 10^^ W Hz“^ and Lr/o/9-4 X 10^^ 
W. 






